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1. Introduction 

Open systems typically give rise to resonances. A resonance is a long-living quasi- 
stationary state, which eventually decays into the continuum. Physically, it may be thought 
of as a particle, initially trapped inside the system, which eventually escapes to infinity. 

One common approach to studying resonances is based on the analytic properties of 
the scattering matrix S(£) in the complex energy plane. Resonances correspond to poles 
£ n = E n — -|r n of S(£) on the non-physical sheet [1,2]. In an alternative equivalent 
approach, one solves the Schrodinger equation subjected to the boundary condition of 
purely outgoing wave outside the range of the potential. This boundary condition, which 
describes a process in which a particle is ejected from the system, renders the problem non- 
Hermitian. The Schrodinger equation with these boundary conditions leads to complex 
eigenvalues £ n which correspond to resonances [1,2]. For a recent lucid discussion of 
resonances in quantum systems, with particular emphasis on the latter approach, see [3,4]. 

The outgoing-wave approach leads, in a natural way, to non-Hermitian effective hamil- 
tonians ,whose complex eigenvalues are the resonances of the studied system [5-7]. Such 
effective hamiltonians are very useful for studying resonances in scattering theory, includ- 
ing scattering in chaotic and disordered systems [8-11]. 

There are many examples of resonances in atomic and nuclear physics. Recently, there 
has been considerable interest in resonances which arise in chaotic and disordered systems. 
See [9] for a recent review. One of the main goals in these studies is computation of the 
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distribution P(T) of resonance widths. There is ample amount of work on computing 
P(r) in one-dimensional disordered chains [10-16]. Numerical results presented in some 
of these works indicate that P(T) ~ T~ 7 in a large range of values of T, where the 
exponent 7 is very close to 1. 

The present work was motivated by [10,11]. In particular, an analytical approach was 
developed in [11] for studying resonances, which is based on counting poles of the resol- 
vent of the non-Hermitian tight-binding effective hamiltonian of the open chain. In the case 
of a semi-infinite disordered chain, coupled to a semi-infinite perfect lead, these authors 
have derived an exact integral representation for the density of resonances (DOR), valid 
for arbitrary disorder and chain-lead coupling strength. In the limit of weak chain-lead 
coupling (in which resonances are typically narrow) they were able to rigorously derive a 
universal scaling formula for the DOR, valid for any degree of disorder and everywhere 
inside the unperturbed energy band of the closed chain. The 1/T behavior of the DOR 
follows from that formula. 

In this paper we shall lay out an analytical approach for computing the DOR of a semi- 
infinite disordered system in the continuum. Our method is based on counting zeros of the 
spectral determinant of the continuous effective non-Hermitian hamiltonian. 

2. Resonances in the Semi-Infinite Disordered System and the Resonance Spectral 
Condition 

Let us consider a disordered system which occupies the domain < x < L. The system is 
assumed to be closed at x = L (eventually, we shall take the limit L — > 00). It is described 
by the Schrodinger hamiltonian 1 

^system plus contact = ~ + ^ ( X ) ~ ^ ' W 

where V(x) is a white-noise potential, drawn from the probability distribution 




Here Z is a normalization constant and D is the variance, i.e., (V(x)V(y)) = D5(x — 
y). The singular last term in (1) is the contact potential, describing the coupling of the 
disordered system to the perfect lead 2 . The lead is described, of course, by the free particle 
Schrodinger hamiltonian i?i ca d = —jd^ along the negative half-line — 00 < x < 0. By 
varying A we can control the coupling of the disordered system to the lead. In particular, 



1 We use units in which m = % = 1. 

2 The hamiltonian (1) can be obtained from the tight-binding hamiltonian used in [10,11] in the 
usual manner, namely, by sending the bulk hopping amplitude t to infinity and the lattice spacing 
a to zero, such that ta 2 — ^— is held fixed. In particular, the contact coupling A is obtained by 
demanding that y = exp 4p ~ 1, where i! is the hopping amplitude associated with the link which 
connects the disordered system and the lead. 
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the limit A — > — oo corresponds to closing the system and disconnecting it from the lead. 
Thus, the DOR should collapse onto the real energy axis and coincide with the density of 
states of the closed disordered system in this limit. 

Resonances in our open system correspond to quasi-stationary solutions of the 
Schrodinger equation 

H^(x,k) = ^(x,k) (3) 

subjected to the outgoing wave boundary condition 

ip(x, k) = e~ ikx , Rcfc>0,lmfc<0 (4) 

in the lead x < [1^,10,1 1]. The condition Rc k > in (4) is imposed because the out- 
going wave propagates to the left in the lead, and Im k < is imposed since the modulus 
of the resonance wave function has to grow into the lead. (Anti-resonances in our system, 
on the other hand, correspond to values of k lying in the third quadrant of the complex k 
plane.) [3,4]. 

In addition to resonances, the open system has also true bound states at negative energies, 
which cannot leak into the lead. They correspond to solutions of (3) subjected to the 
outgoing wave boundary condition (4) in the lead, at values k which lie on the positive 
imaginary axis (their wave functions therefore decay into the lead). 

Motivated by [17,18], it is useful at this point to introduce the logarithmic derivative 

'<*•*> -^f- 

It follows from (1) and (3) that inside the system f(x,k) obeys the Riccati equation 

d x f + f = 2V(x) -k 2 , 0<x<L, (6) 

where f(x, k) is subjected to the initial condition f(L, k) = — oo, independently of k. This 
is due to the Dirichlet boundary condition 3 set upon ip(L, k), since the system is closed at 
x = L. In practice, we shall impose the initial condition 

f(L,k) = -c (7) 

for some real constant c, and in the end take the limit c — > oo. 

Evidently, in the lead, f(x, k) is constant for any solution of (3) that is subjected to the 
outgoing wave boundary condition (4). In particular, for such a solution, 

f(0-,k) = -ik. (8) 

Therefore, due to the 5(x) term in (1), we know that the logarithmic derivative has to jump 
to 



3 ip{x,k) either vanishes at x = L through positive values and nonvanishing negative slope, or 
through negative values and positive slope. 
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f(0+,k) = -ik-\ (9) 

right across the contact, on the system's side. We see that in the limit A — > — oo, we are 
effectively imposing a Dirichlet boundary condition on tp(0+, k), which means closing 
the disordered system and disconnecting it from the lead, in accordance with the comment 
made above Eq.(3). 

We now have all the ingredients required to solve for resonant states in our system. For 
a given instance of the random potential V(x) drawn from (2), and for a given complex 
wave-number k, we integrate the Riccati equation (6), subjected to the initial condition (7) 
backwards in x, all the way down to x = 0+. Note that due to this integration backwards, 
f(x, k) can only depend on values of V(y) with y > x. More precisely, it is straightfor- 
ward to deduce from (6) that 

F(x, y) = = -2% - x) cxp -2 J V f(z, k) dz . (10) 

(We shall make use of this result in the next section.) Once we have /(0+, k), we compute 

m M) - /(0 + ,t) +ik + \= M»MHt» + W»» , (U) 

whose functional dependence on V{x) is indicated explicitly. Thus, from (9), 

F(k;[V])=0 (12) 

whenever we have a resonance at k (in which case k lies in the fourth quadrant, off the 
coordinate axes), an anti-resonance (k e third quadrant, off axes), or a bound state at 
negative energy (in which case k = i\k\ is pure positive imaginary). 

Note that F(k; [V]) has also poles, which occur whenever 4 ?/' (0, k) — 0, i.e., whenever 
tp(x, k) obeys the Dirichlet boundary condition at x = 0, in addition to the boundary con- 
dition d x ip(L, k) + cip(L, k) = it is already subjected to from (7). Thus, these zeros 
correspond to eigenstates of the Hermitian system (closed at x = 0), whose energy eigen- 
values E — i/c 2 can be either positive (with k real), or negative (with k pure imaginary, 
on the positive imaginary axis). It is the positive-energy eigenstates of the Hermitian dis- 
ordered system which become resonances upon connecting it to the lead at x = 0. Its 
negative energy eigenstates, while suffering slight distortions and energy shifts due to the 
contact potential, cannot leak out into the lead, and remain bound states of the open system, 
as was mentioned above. 

To summarize, poles of F(k; [V]) lie exclusively on the real and imaginary axes of the 
complex k -plane. They correspond to eigenstates of the closed Hermitian system. Zeros 
of F(k; [V]) correspond either to resonances, anti-resonances, or to genuine negative en- 
ergy bound states of the open system. We thus conclude that F(k; [V]) is essentially the 
ratio of two spectral determinants: the spectral determinant of the non-Hermitian hamilto- 
nian, which gives rise to resonances, divided by the spectral determinant of the Hermitian 
problem [19]. 



4 As indicated explicitly in (11), F{k; [V]) is the ratio of two functions. These are analytic functions 
of k. Therefore, poles can arise only as zeros of the denominator. 
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The quantity F(k; [V]) is a holomorphic function of k, which follows immediately from 
the holomorphy of wave function tp(x, k). Counting resonances in our system is thus 
equivalent to counting the zeros of F(k; [V]) in the appropriate region in the complex k- 
plane, namely, the fourth quadrant (and off the coordinate axes). To this end, we follow 
[20,21]. Let ko be one of the zeros of F(k; [V]). For simplicity, we shall take it to be a 
simple zero 5 . Thus, in the vicinity of ko, F(k; [V]) ~ F'(ko)(k — ko), and therefore 6 

^a& log|F(fc;[t/])|2 = 5(2)(fc - fco) ' (13) 

where we used the identity d ^ dk log \k\ 2 = nS^ (k) . Similarly, a pole in F(k; [V]) at 
k = ko (which must lie either on the real or the imaginary axis), would result in (13) with 
a negative sign on the RHS. 

Off the coordinate axes, F(k; [V]) has only zeros, not poles. 
Thus, dk*dk log \F(k; [V])\ is a sum of strictly positive terms of the form (13). Let us 
denote the set of all these zeros of (11) by S. Therefore, we obtain the density of these 
zeros as 

p(k, k*, [V]) = 5 ( Rc ( fc - M) 8 (Im(fc - k n )) 

where k is kept off the coordinate axes. In the fourth quadrant Refc > 0, Imk < 0, the 
distribution p(k, k* , [V]) coincides with the DOR of our system, for the given instance 
V(x) of the disorder potential. 



3. The Average Density of Resonances 

We would like now to average (14) over the disorder (2), in order to obtain the average 
DOR 

p(k,k*) = I^^log|F(fc;[y])| 2 ), Refc>0,lmfc<0 (15) 

of the disordered system. To this end we need the probability distribution function (pdf) of 
the random variable /(0+, k), which enters (11). In other words, we need the pdf for the 
event that /(0+, k) = u + iv, namely, 

P(u, v) = (Re/(0+, k)-u)5 (Im/(0+, k) - v)^. (16) 



5 We lose nothing by restricting to simple zeros. Indeed, if fco is a multiple zero of order n, then 

F(k; [V]) ~ F J, k °\ k — fco)™, and as should be clear from the discussion below, it will be simply 
counted n times. 
6 Here <5 (2) (fc) = c5(Refc) 6 (Imk). 
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(In order to avoid notational cluttering, we have suppressed the explicit dependence of this 
distribution on k,k*.) Then, from (1 1), (15) and (16), we obtain the integral representation 



1 d 2 



u + iv + ik + X 



2 



(17) 



for the average DOR, in terms of the pdf P(u, v) for /(0+, k). Eq. (17) is our main result. 

Thus, we need to determine P(u, v). This we do by extending the well-knwon treatment 
of the disordered closed system in [17] (which is neatly explained and summarized in [18]), 
to the case of open systems. Thus, we introduce the pdf 

P(u, v; x) = (5 {Ref(x, k) - u) 5 (Im/(a;, k)-v)} (18) 

for the event that f(x, k) = u + iv. We observe that the Riccati equation (6) for f(x, k) can 
be interpreted as a Langevin process in the complex plane, with the coordinate x playing 
the role of "time". By standard methods of the theory of stochastic processes (which we 
do not belabor here), we can derive the Fokker-Planck equation 

d x P = -2Dd\P + d u [{u 2 - v 2 + 2E r )P] + d v [(2uv + 2Ei)P] (19) 

for P(u, v; x). Here E r and Ef are, respectively, the real and imaginary parts of the com- 

,2 

plex energy E = 4p We should perhaps just mention that in deriving (19) we used the 
relations 

= -20(0) = -1 and Im^ = , 

which follow from (10). 

A couple of remarks on (19) are in order at this point: First of all, note the diffusion 
term in (19). The minus sign in front of d 2 P arises because "time" runs backwards. The 
diffusion constant is of course positive. Second, note that there is no diffusion term in the 
v direction. The reason for this is that V(x) is strictly real. If we allow also a white noise 
imaginary part for V(x) with variance e, which is uncorrected with its real part, we will 
have a diffusion term —2ed 2 in the v direction as well. In fact, it can be shown that the re- 
sulting Fokker-Planck equation in this case is related, via a similarity transformation (i.e., 
by a gauge transformation with an imaginary vector potential) to a Schrodinger equation, 
in euclidean time, for a particle with anisotropic mass (assuming D ^ e), with an inhomo- 
geneous magnetic field B oc v perpendicular to the uv plane, and a quartic potential. 

Let us now resume our discussion of P(u, v). By integrating (19) backwards, from the 
initial point x = L, where (7) implies the (k independent) initial condition 

P(u, v- L) = 5(u + c) 5{v) , (20) 

all the way down to x = 0, we can determine the desired pdf P(u, v; 0) = P(u, v). In the 
limit of very long system L — > oo, the pdf at x = 0, P(u, v), which goes into (17), should 
be governed by the simpler stationary Fokker-Planck equation 

-2Dd 2 u P + d u [(u 2 -v 2 + 2E r )P] + d v [(2uv + 2Ei)P] = . (21) 
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It is straightforward to check, directly from (21), that the leading asymptotic behavior of 

P(u, v) as r 2 = u 2 + v 2 — > oo is 



(up to overall normalization). Thus, P(u,v) decays faster than the minimal rate required by 
normalizability, and in a manner which is independent of k. In fact, more careful analysis 
shows that k (or energy) dependence appears for the first time only in the r -6 term. 

Note that in the case of a closed disordered system, discussed in [17,18], with Dirichlet 
boundary conditions at both x — and x = L, f(x, k) is purely real (all eigenstates 
tp{x, k) are real, up to an overall phase). Thus, all reference to v, the value of Im/(;r, k), 
disappears from (19) and (21). In particular, the stationary Fokker-Planck equation in this 
case is simply 



with real energy E, and which can be solved explicitly, of course, by elementary methods. 
At the moment, we are making an effort to derive a certain integral representation for 
P(u, v), the solution of (21), which we shall report on in a subsequent publication. 
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